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Abstract. We present results for a finite variant of the one-dimensional Toom model with 
closed boundaries. We show that the steady state distribution is not of product form, but is 
nonetheless simple. In particular, we give explicit formulas for the densities and some nearest 
neighbour correlation functions. We also give exact results for eigenvalues and multiplicities of 
the transition matrix using the theory of ^-trivial monoids in joint work with A. Schilling, B. 
Steinberg and N. M. Thiery. 


1. Introduction 

Dynamical variants of the Ising model have been of interest since Glauber’s pioneering work [T] 
on the subject. Although there are many dynamical rules for the Ising model whose steady state 
distribution is Gibbsian, there are many other dynamical rules whose steady state distribution 
has very different structure, see for example 

The model that we will be interested in is a probabilistic cellular automaton of Ising spins 
on Z 2 due to A. Toom |4J. The dynamics is given by the so-called NEC (north-east-center) rule. 
Recall that cellular automata are discrete-time dynamical models, where the entire configuration 
is updated simultaneously. Denote the spin at the site (i,j) at integer time t E Z + as <Tjj(t), 
where aij(t) can be ±. The dynamics then evolves as follows. 

+ with prob. p 

— with prob. q 

sign(crjj(t) + + a itj+1 (t)) with prob. 1-p-q 

To be more explicit, the spin at a given site at time t + 1 evolves according to three spins at 
time t; itself, its north neighbour and its east neighbour. We illustrate the rule below, where we 
only show the case where cnj(t) = +. The other case is similar. For clarity, we have circled the 
changing spin. 


o i,j{t + 1) = < 


+ v + + o + — l-p — 

©+ ©+’ ©- ©-’ ©- © - 


Notice that the model becomes a deterministic NEC majority rule when p = q = 0. lip ^ q, the 
rule favours one spin over the other, which can be thought of as the effect of a magnetic field. 




It turns out that the model exhibits a noise-dependent phase transition, which was 
demonstrated rigorously by Toom [3], He showed that for sufficiently high p, q, there is a 
unique steady state, whereas for low enough p, q, there are at least two steady states. 

Derrida, Lebowitz, Speer and Spohn 00 considered this model on the third quadrant in Z 2 , 
with the boundary condition that all spins on the negative x-axis are — and all on the negative 
y-axis are +. When p = q = 0, any configuration in which +’s and — ’s are separated by a single 
interface of staircase shape is stationary, see Figure |T] They were interested in the dynamics 
of the interface for very small p, q. First, notice that in that case, all defects in the bulk of 
these phases are quickly cancelled out. However, the interface can and does fluctuate because 
of spontaneous sign changes of the +’s below and the —’s to the left of the interface. Figure Q] 
shows how this happens for a + spin below the interface. Notice that the new interface is still 
of staircase shape. 
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Figure 1. Dynamics in the PCA at low temperature in the third quadrant. The + spin which 
will spontaneously flip is shown on the left as © and the spins which flipped to — are shown in 
the right as ©. 


Derrida, Lebowitz, Speer and Spohn 00 understood the interface dynamics by thinking of 
the interface as itself an Ising spin configuration S on Z + as follows. Look at the steps taken 
by the interface starting from the north-east corner. If the n’th step is vertical, set S n = +, 
otherwise set S n = —. 

The transition shown in Figure [T| corresponds to the transition 

4— -1— I-1—I—I—h • • • —4— 4 - 1-1—I—I—b • • • , (LI) 

where the part undergoing the change has been underlined. In this one-dimensional model, 
each ± exchanges with the first © on its right with rate X±. They also proposed to call this 
one-dimensional version the “Toom model”. 

However, there turns out to be a simplification of this model due to Lebowitz, Neuhauser 
and Ravishankar (LNR) [7j. Rather than allowing all spins to exchange, they only allow the 
leftmost spin in a block to exchange with the first opposite spin to its right. The transition in 
m would not be allowed in this model. An example of an allowed transition is 


+-+ + - + + + + ••• 


++-+ - + + + + ■■■. 


( 1 . 2 ) 






















Further, the first spin flips independently with a prescribed rate a. They studied the steady 
state distribution of this model when A_|_ = A_ = a = 1. The most interesting result there from 
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our point of view is [?., Theorem 2] that the density at site n goes like __ 

It is this variant that we will refine. We consider the model on an interval of L sites with 
closed boundaries. We modify the LNR model by introducing two parameters a and f3. In the 
original model [?), a and (3 were taken to be 1. We will also make the connection to exclusion 
processes explicit by replacing the spins ± by 0,1. Each site is then occupied either by a particle 
of type 0 or that of type 1. The dynamics is as follows: the leftmost particle in a block of particles 
of the same type can exchange with the opposite type particle to the right of the block. Clearly, 
the site preceding the block must be of opposite character if there is such a site. This exchange 
occurs with rate a if the block contains 0’s and /3 if the block contains l’s. Note that the block 
size can be one. To recap, 


0...01 



k 




with rate a, 
with rate (3, 
with rate a, 
with rate (3 , 


if the block begins at the first site, and 


if the block does not begin at the first site. 


(1.3) 


An example of this dynamics is given with different notation in (11.21) . The model is particle 
conserving and one defines it using integers (no,ni), with no + n\ = L, to denote the number 
of holes and particles respectively. There are thus ( n °,+ ni ) configurations. We will call this 
the (no, ni)-system. The rules then are just given by (11.31) for any positive value of k. In 
particular, the rightmost block cannot by itself make any transitions. It is of course affected by 
the transitions involving the block immediately to its left. 

Note that the presence of the boundary is crucial to having a nonequilibrium steady state. If 
we consider the model on a ring of L sites, then there is a unique absorbing set of configurations 
with all no holes in one block and all n± particles in another block. These configurations are 
translations of one another, and the steady state is uniformly distributed on this small set of 
configurations. On the interval model, however, the boundary acts to break up the size of blocks. 
For example, in the (3,3) system, if one starts with the configuration 001110, then one can have 
a transition to the configuration 100110 with rate a. Then one has split the l’s from a block of 
size 3 to two blocks. 

We will present results for this model in the rest of the paper. The proofs will be delegated to 
a subsequent publication [8] . In Section O we demonstrate the structure of the Markov matrices 
for the model. We will first give a recurrence for constructing the matrices, and give formulas for 
their eigenvalues. In Section [3l we present formulas for the density and other nearest neighbour 
correlations. We will also present a conjecture for the partition function. 


2. The Markov Matrices 

Here we will show that the Markov matrices for the model have a rich structure. Recall that 
the Markov matrix of the (no, ni)-system is a matrix of size ( n ° r ~ ni ) indexed by configurations 
of the system whose (i,j)’th entry is the rate of going from state j to state i. We will denote 
the matrix by M n0)ni . 

As an illustration, consider the example with no = ni = 2. The configurations, in 
lexicographic order are: (0011,0101, 0110,1001,1010,1100). The Markov matrix in this ordered 



basis is given by 


M 2,2 


/ —a /3 P 

0 -p-2a 0 

0 a — P — a 

a a 0 

0 0 a 

0 0 0 


0 0 0 \ 

p 0 0 

0 P P 

~P — a P 0 

0 -2 p -a 0 

a a ~P ) 


( 2 . 1 ) 


The steady state distribution is the column (or right) eigenvector with eigenvalue 0. We will 
normalise it so that the sum of the entries is 1, i.e. it represents a probability distribution. We 
will denote the vector by Tr no , ni - For this case, 

P 2 a.p 3 a 2 P (a + 2P) aP 2 (2 a + P) a 3 P a 2 

(a + P) 2 ’ (a + P) 4 ’ (a + /3) 4 ’ (a + /3) 4 ’ (a + /3) 4 ’ (/3 + a) 2 

Note the simple probabilities of configurations with just two clusters. We will also see later that 
M no m has simple formulas for eigenvalues and their multiplicities. The eigenvalues of M 22 are 




0 with multiplicity 1, 

—a — P with multiplicity 4, (2-3) 

—2a — 2 P with multiplicity 1. 


We will begin by explaining the recursive structure of the Markov matrices. To that end, we 
will first need to define some auxiliary matrices. Let 1*. denote the k x k identity matrix and let 
A no , ni be a matrix of size x ( n °" l ^ L 1_1 ) defined as a concatenation of identity matrices 

of increasing sizes which are justified at the top row as, 


*-no,ni 


l/ni-n 1/ n ± \ 


n -0+ n l~ 2 \ 

ni~l ) 


(2.4) 


and the rest of whose entries are zero. For example, 



/ 1 

1 

0 

0 

1 

0 

0 

0 

0 

0 ^ 



0 

0 

1 

0 

0 

1 

0 

0 

0 

0 



0 

0 

0 

1 

0 

0 

1 

0 

0 

0 


^3,3 = 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

1 

0 




0 

0 

0 

0 

0 

0 

0 

0 

1 ) 


a matrix of 

size 

/no+ni 
V n 0 — 

1 J 

1 X 

/no+ni 
V n 0 

-1 

) defined 


Similarly, let B n0)Ul 
identity matrices of increasing sizes which are justified at the bottom row, 


Bn 0 ,n 1 — ( 

1/n 

V ^ 

o+™ 
n o - 

l-2> 
-1 ) 

1 


1 

( n 0 
\n Q - 

1 ) 

1 

1 1 

0 0 
e_5 

and the rest of whose entries are 

zero. 

For example, 






( 1 

0 

0 

0 

0 

0 

0 

0 

0 

0 ^ 



0 

1 

0 

0 

0 

0 

0 

0 

0 

0 


F>3,3 = 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

1 

0 

0 

1 

0 

0 

0 



0 

0 

0 

0 

1 

0 

0 

1 

0 

0 



v° 

0 

0 

0 

0 

1 

0 

0 

1 

1 ) 





We will write down the Markov matrix in the lexicographically ordered basis of configurations, 
just as in the example (12.111 . The transition matrix M no ^ ni can be expressed in terms of M no _i ni 
and M n0jni _i using 2x2 blocks as, 


M, 


no,n i 




IWM’O) 

M(M) 


where each block can be expressed in terms of smaller blocks as 


M(°-°) = M no _ 1>ni - a (■ 
M (oa) = 0 


A, 


no,ni 


M (h°) = 


a 


no,ni 


no,ni 


M&V = M no , ni _i - P ( Bno 0 ni ) 


Here 0 represents a zero matrix, A nQyU1 and B n0in 1 are defined in ((2.4(1 and (12.5(1 respectively and 
with the initial condition Mq u1 = M n0i o = (0), the lxl zero matrix. The sizes of the blocks 
of zero matrices is fixed by the fact that M ( 0, °) and IW^ 1 ’ 1 ) are square matrices of size ( n °^ ni ) 
and ( n °+ n 1 ) respectively. One can check that the matrix of the example at the beginning of this 
section in (12.1(1 can be recovered from the recursion above. The proof of this recursion follows 
from an elementary analysis of transitions in the ordered basis. For instance, M ( 0,1 1 represents 
all transitions from configurations that begin with a 1 to those that begin with a 0. 

The deepest results that have been obtained about this model are on the eigenvalues and 
their multiplicities of the Markov matrices in joint work with A. Schilling, B. Steinberg and N. 
M. Thiery [9]. The result there is more general, but it reduces to the following in our case. 
The characteristic polynomial of the transition matrix M no>ni is given by the explicit product 
formula 

min(no, ni) 

(n 0 \ 

|A1 — M no>ni | = (A + k(a + P)) \ k A k) , (2.6) 

k =o 


where the factor for k = 0 corresponds to the zero eigenvalue, which naturally occurs with 
multiplicity one. One can check that the eigenvalues for no = n\ = 2 in (|2.6I) match those given 
in (12.31) . 

The proof of (12.61) is highly nontrivial and has been obtained using techniques from the 
representation theory of ^-trivial monoids; another application of this theory is to a nonabelian 
model of sandpiles m- No elementary proof of ((2.61) has been found so far. 


3. Correlation Functions 

We will now state exact results for some correlation functions including the density. In some 
cases, the analysis is straightforward and will be explained, but the more complicated cases will 
be stated without proof. Detailed proofs will appear elsewhere [8]. 

We will use r}i to denote the occupation variable of the particle at site i and angle brackets (•) 
to denote averages in the steady state distribution. As before, we will set uq + n\ = L. First of 
all, note that the system satisfies an obvious particle-hole symmetry. Namely (?yi,..., rn) in the 
(no, ni)-system is equal to (1 — 771 ,..., 1 — tjl) in the (n\, no)-system with a and P interchanged. 










It turns out that the nature of explicit formulas for correlation functions is different for the 
sites less than both no and ni. Consequently, it will be useful to denote n* = min(no,ni). This 
is because looking at the conditional process on the first n* sites is equivalent to looking at the 
first n* sites in the semi-infinite system [7] because no matter how large a block of particles 
(holes) is within this subsystem, there is always a hole (particle) far enough to the right. 

First of all, the density at the first site is given by 


(m) 


a 

ol + f3 


(3.1) 


This is clear since there is a hole somewhere in the system for an exchange with the leftmost 
site (which is of course the leftmost site for the first block), 




With just a little more work, we can also show that the density in the first n* sites is constant, 
i.e., 


(Vk) 


a 

a + (3 ’ 


(3.2) 


for 1 < k < n*. To see this, suppose that 2 < k < n*, and consider the master equation for (vk)- 
First, we look at the outgoing contribution. If r/ is a configuration with % = 1 then, if r/k-i = 1 
then 7]^ cannot change instantaneously. If r/k-i = 0, then % can become 0 in two ways: at rate 
a due to an exchange with Vk-j f° r some j > 1 and at rate (3 due to an exchange with for 
some j > 1. Thus the total negative contribution to ^{Vk) is 


-(« + /3)(( 1 ~m-i)Vk)- 


Similarly, if i] is a configuration with % = 0 then, if %_i = 0 then % cannot change 
instantaneously. If %_i = 1 then % can become 1 in two ways: at rate (3 due to an exchange 
with r/k-j for some j > 1 and at rate a. due to an exchange with r/k+j for some j > 1. Thus the 
total positive contribution to ^(Vk) is 

(« + /?)(%-i(l - Vk))- 


So 


= (a + P)[-{(l-rik-i)Vk) + {vk-i{3 ~ Vk))} 
= (<x + P)[(Vk- 1) - (Vk)], 


(3.3) 


and this must vanish in the steady state. Using the initial condition in (13.11) proves (J3.2D . 

The reason the above proof does not go through if k > n* is that it assumes there is a particle 
of the opposite type far enough to the right of the block starting at k. When k is too large, such 
an assumption does not hold. The density at a general site is given by the more complicated 
formula, 


(Vk) 


/3 k ~ ni 
(a + (3) k 



(3.4) 


A little bit of manipulation shows that (13.4ji reduces to (|3.2I) when k < n*. The proof of this 
result is more complicated and uses results for other correlations stated in ()3.6j) and (|3.7j) . 






For more complicated correlation functions such as block size distributions, we have some 
preliminary results. For k <n\, 


To see this, first consider the outgoing transitions for a configuration 7] with rji = ■ ■ ■ = rj^ = 1. 
Note that 7 / 2 ,,r)k cannot change instantaneously. The only position which can effect a 
transition is the first one, where a hole to the right of k exchanges with the particle at the 
first site with rate fd. The only incoming transitions to T] are from configurations rf in which 
?/] = ■ ■ ■ = r}' k _i = 1 and r/ k = 0 is nonzero. There is a transition to ij at rate a which exchanges 
a particle to the right of k with the hole at the kth site. Therefore 


(3.5) 


0 = 


d_ 

dt 


id l • • • Vk) 


a(m ■ ■ ■ Vk- 1(! - Vk)) - P(Vi ■ ■ ■ Vk), 


which implies the recursion, 


(m ■ ■ ■ Vk) 


a(Vi ... 7? fc _i) 
OL + fd 


Using (13.11) as the initial condition gives the desired result. 

Our final results are about correlations at a block of sites at the end of the system. The 
probability of having a block of L + 1 — k l’s starting at position k is given by 


{Vk ■ ■ ■ Vl) 



(3.6) 


The formula becomes much simpler if we also look at the joint correlation with a hole at site 


k- 1. 


((1 - Vk-i)Vk ■ ■ ■ Vl) 


( k-2\a k - 1 - n °/3 n ° 
\n 0 - 1/ (a + fd) k ~ l ' 


(3.7) 


Note that both formulas (13.61) and (13.71) give 0 when we set k = no as expected. These will be 
proved in a subsequent publication [8]. One can get the joint correlation of blocks of holes in 
the first and the last sites by interchanging a and fd in (13.51) . (13.61) and (13.71) . We do not yet 
have exact formulas for block correlations in the bulk. 

Steady state probabilities are rational functions of a and fd , and hence there is no unique way 
to write them. Suppose we write the steady state probabilities such that the greatest common 
factor of the numerators is 1. This forces them to be written in a unique way. See for example, 
m - We will call the normalization factor, denoted Z(no,ni), in the (no, ni)-system to be 
the least common multiple of the denominators of the steady state probabilities when they are 
written in this way. The formula for the normalization factor Z(n o, n\) seems to be tantalisingly 
clean, 

Z(n 0 ,n 1 ) = {a + P) noni . (3.8) 


When ?io = n\ = 2, this matches with ()2.2I) . Although this formula seems quite simple, the 
proof has resisted the best of our efforts, and we leave it as a conjecture for now. 

To conclude, let us compare some of our results with those of Lebowitz, Neuhauser and 
Ravishankar |7j. First of all, our model is on the finite one-dimensional lattice, whereas theirs is 
on Z>o- We have a more general model in the sense that we have more parameters. Note that 
the ratio of a and fd is the only extra parameter here, but in our generalised model [9], we have 
many more. We have explicit formulas for many correlation functions in (I3.4|) . (l3.5p . (l3.7p . (l3.6p . 
and a conjecture for the partition function (I3.8p . 







One can compare the densities in both models by performing asymptotics of (13.41) . Of course, 
when a = (3 and k is large, the density at site k goes exactly like —exactly as in [7j. A 

Z\ 7TrC 

preliminary analysis seems to indicate that the density decays like 1 jk for all values of a//3. 
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